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ABSTRACT 


A  study  was  conducted  to  develop  a  mathematical  model  of  a  jet 
exhausting  from  a  flat  plate  into  a  crossflowing  stream.  The  modeling 
was  accomplished  using  the  vortex-lattice  method.  Analytical  stream¬ 
lines  and  pressure  distributions  on  the  flat- plate  surface  as  well  as 
vector  data  in  the  field  above  the  plate  are  compared  with  available 
experimental  data.  Recommendations  are  made  for  further  improve¬ 
ment  of  vortex- lattice  jet- modeling  techniques. 
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NOMENCLATURE 

Cp  Pressure  coefficient,  (p_pa>)lq 

D  Jet  exit  diameter,  ft 

P  Pressure,  lb/ f t^ 

o 

q  Dynamic  pressure,  lb/ ft^ 

R  Jet  exit  radius,  ft 

V  Velocity,  ft/ sec 

Ve  Effective  velocity  ratio, 

X,  Y,  Z  Coordinates  related  to  free- stream  velocity 

6  Deflection  of  jet- exit  centerline,  positive  with  increasing 

angle  away  from  free- stream,  deg 

p  Density,  slugs/ft^ 

SUBSCRIPTS 

j  Jet  exit 

®  Free  stream 
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SECTION  I 
INTRODUCTION 


Much  work  has  been  done  to  define  the  characteristics  of  jets  ex¬ 
hausting  into  a  crossflowing  stream.  Numerous  experimental  programs 
have  been  performed  to  gain  insight  into  the  behavior  of  these  jets  and 
their  mutual  interaction  with  the  surrounding  fluid.  Schemes  for  ana¬ 
lytically  duplicating  the  complex  flow  fields  associated  with  jets  have 
been  developed,  with  varying  degrees  of  success,  by  a  number  of  indi¬ 
viduals. 

This  report  presents  the  results  of  an  attempt  to  apply  a  previously 
developed  theory  for  representing  jet  plumes  to  a  configuration  in  which 
the  jet  exhausts  from  a  flat  plate.  The  technique  presented  allows  com¬ 
putation  of  pressures  and  streamlines  on  the  flat  plate.  Velocity  vectors, 
streamlines,  and  pressures  can  be  calculated  throughout  the  flow  field 
surrounding  the  jet  plume.  The  method  does  not  require  a  great  deal  of 
computer  time,  even  though  it  is  capable  of  handling  complex  configu¬ 
rations  of  multiple  jets  and  solid  surfaces. 


SECTION  II 

PROPERTIES  OF  THE  JET 


The  vortex- lattice  model  described  in  this  paper  is  an  attempt  to 
simulate  the  physical  phenomena  observed  in  jets  exhausting  into  cross- 
flows.  Consideration  is  given  to  representing  the  effect  of  the  jet  upon 
the  surrounding  flow  field  and  the  flat-plate  surface.  Simulation  of  the 
internal  flow  is  ignored.  Therefore,  it  is  necessary  to  have  a  basic 
understanding  of  the  properties  of  such  jets  and  the  associated  aero¬ 
dynamic  characteristics. 

The  jet  exiting  a  surface  enters  the  stream  normally  (i.  e. ,  without 
the  skew  angle  exhibited  in  the  free- jet  case).  As  the  jet  exhausts  from 
the  surface,  it  is  bent  in  the  direction  of  the  crossflowing  stream  be¬ 
cause  of  two  primary  mechanisms.  The  major  cause  is  the  momentum 
transfer  due  to  turbulent  shear  forces  acting  on  the  lateral  regions  of 
the  jet.  This  results  in  general  sweeping  downstream  as  fluid  from 
the  crossflow  is  entrained  into  the  jet.  A  smaller  influence  on  the  jet 
bending  results  from  the  pressure  differential  across  the  jet.  Several 
investigations  have  been  conducted  to  define  the  jet  deflection,  and 
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numerous  empirical  relationships  have  been  derived  for  the  axial  cen¬ 
terline  trajectory.  All  agree  that  the  jet  path  trajectory  is  mainly  a 
function  of  the  velocity  ratio.  The  expression  used  in  this  development 
was  provided  by  Margason  (Ref.  1).  as  follows: 


7_ 
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Since  the  sides  of  the  jet  are  deflected  downstream  more  rapidly 
than  the  core,  the  cross  section  takes  on  the  classical  crescent  shape. 
This  eventually  results  in  a  rollup  of  the  jet  and  formation  of  two  contra¬ 
rotating  vortices  which  endure  relatively  far  downstream  and  then  fi¬ 
nally  submit  to  viscous  effects. 

As  the  jet  penetrates  the  crossflowing  stream,  its  cross  section 
increases  in  size  because  of  both  expansion  and  the  addition  of  free- 
stream  fluid.  The  two  mechanisms  responsible  for  this  addition  of 
fluid  have  been  dubbed  turbulent  and  nonturbulent  entrainment.  Turbu¬ 
lent  entrainment  takes  place  in  the  regions  of  turbulent  mixing.  This 
is  the  same  means  by  which  jets  exhausting  into  quiescent  air  entrain 
and  is  caused  by  mean  velocity  differences  within  the  flow.  Nonturbu¬ 
lent  entrainment  is  primarily  a  result  of  the  contrarotating  vortex  pair. 
Free- stream  fluid  is  swept  around  the  sides  of  the  jet  and  is  captured 
by  the  jet  rollup  process.  Previous  studies  indicate  that  entrainment 
caused  by  the  latter  pressure- induced  phenomena  may  be  an  order  of 
magnitude  greater  than  that  caused  by  the  viscous  mechanism. 

A  jet  exhausting  into  a  crossflow  displays  a  low-pressure  wake 
region  on  the  downstream  side.  It  has  been  shown  that  the  jet  exhibits 
a  blockage  similar  to  a  cylinder  with  suction.  As  free-stream  fluid 
approaches  the  jet,  it  is  first  influenced  by  the  blockage  effect.  It  is 
deflected  outward  and  upward  by  the  contrarotating  vortices  and  is 
either  captured  by  one  of  the  entrainment  mechanisms  or  swept  into  the 
aft  wake  region.  A  number  of  plots  of  the  q-vector  data  given  by  D.  K. 
Mosher  (Ref.  2)  are  shown  in  Fig.  1,  Appendix  I.  Note  that  without  a 
knowledge  of  the  total  pressure  values  it  is  impossible  to  distinguish 
the  jet  from  the  free  stream,  since  sharp  discontinuities  in  flow  direc¬ 
tion  occur  only  near  the  jet  exit. 

A  major  portion  of  jet  investigations  have  dealt  with  the  pressure 
field  induced  on  the  surface  from  which  the  jet  exhausts.  In  V/STOL 
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applications,  a  primary  concern  is  the  lift  loss  caused  by  the  jet.  Prac¬ 
tical  application  of  an  analytical  model  requires  accurate  simulation  of 
this  phenomenon. 

Pressures  on  the  surface  are  known  to  be  affected  by  a  combination 
of  blockage  and  entrainment  mechanisms.  Mosher's  study  indicates  that 
both  mechanisms  exert  influence  at  low  velocity  ratios;  however,  en¬ 
trainment  effects  dominate  as  the  velocity  ratio  increases.  Blockage 
affects  the  areas  immediately  to  the  front  and  rear  of  the  jet,  whereas 
entrainment  effects  encompass  all  the  remaining  areas.  Plate  pres¬ 
sure  contours  from  Ref.  2  are  shown  in  Fig.  2  for  a  variety  of  velocity 
ratios.  As  the  flow  approaches  the  jet,  it  is  decelerated  because  of  the 
blockage,  and  a  small  positive  pressure  area  is  formed.  The  flow  is 
then  diverted  around  the  jet  and  is  accelerated,  resulting  in  a  negative 
pressure  region  on  the  sides  of  the  jet.  The  contrarotating  vortices 
have  a  major  effect  on  the  plate  flow  as  they  entrain  fluid.  They  cause 
an  insweep  of  flow  to  the  sides  of  the  jet,  resulting  in  a  sink-like  pat¬ 
tern  which  causes  negative  pressure  regions  far  from  the  jet  sides. 

When  the  velocity  ratio  is  increased,  the  entrainment  mechanism  be¬ 
comes  stronger.  The  positive  pressure  region  forward  of  the  jet  is 
diminished  in  both  size  and  magnitude,  and  the  negative  areas  to  the  sides 
and  rear  expand  and  become  more  negative.  Finally,  at  some  velocity 
ratio,  the  positive  region  will  disappear.  It  should  be  noted  that  sev¬ 
eral  authors  have  carried  out  similar  investigations  and  have  arrived 
at  the  same  general  conclusion  (Refs.  2,  5,  6,  and  7).  Data  from 
these  experiments  are  available  in  the  literature. 

Additional  information  concerning  the  flow  over  the  plate  surface 
can  be  obtained  through  the  use  of  oil-film  traces.  The  result  of  an 
oil  study  conducted  by  Mosher  appears  in  Fig.  3  along  with  the  authors' 
interpretation  of  the  plate  surface  streamlines.  As  the  figure  shows, 
the  flow  into  the  region  aft  of  the  jet  exit  is  not  continuous.  A  pair  of 
stagnation  lines  is  shown  extending  downstream  from  the  back  of  the 
exit.  It  is  evident  that  the  lines  result  from  a  pair  of  contrarotating 
eddies  which  lie  very  near  the  plate  surface.  The  disturbances  from 
these  eddies  are  not  believed  to  extend  far  up  into  the  wake  region, 
although  no  quantitative  flow-field  data  could  be  found  in  the  literature 
defining  the  volume  in  the  vicinity  of  the  plate.  However,  the  vector 
information  recorded  in  Ref.  2  shows  that  no  disturbance  occurs  be¬ 
yond  two  diameters  above  the  surface.  Since  the  analytical  model 
described  in  this  document  provides  no  means  for  representing  these 
eddies,  the  stagnation  lines  are  not  expected  to  appear  in  the  com¬ 
puted  flow  pattern  on  the  plate. 


3 


A  E  D  C-TR  -73-57 


SECTION  III 

JET  MODEL  CONSTRUCTION 


Reference  3  presents  the  development  of  a  general  jet  efflux  repre¬ 
sentation  by  the  vortex-lattice  method.  This  jet  model  is  composed  of 
an  inlet  simulation  and  exhaust  blockage  represented  by  lifting  vortex 
panels.  Given  the  velocity  ratio,  a  computer  program  automatically 
computes  the  trajectory,  inlet  skew  angle,  and  tube  geometry.  This 
model  was  used  to  simulate  a  free  propeller  submerged  in  a  crossflow. 

Preparation  of  an  analytical  model  for  a  jet  exhausting  from  a  flat 
plate  required  four  basic  changes  to  the  previously  developed  propeller 
jet  representation.  The  vena  contracta  associated  with  the  propeller 
was  removed,  leaving  a  vortex- lattice  tube  of  uniform  cross  section. 

The  tube  axis  was  forced  to  intersect  the  jet-face  plane  at  an  angle  of 
90  deg,  rather  than  at  the  skew  angle  computed  from  the  velocity  ratio 
(i.  e. ,  Vj /V^).  In  the  initial  configuration,  a  large,  flat  vortex  grid 
was  attached  to  the  jet  face  to  represent  an  infinite  flat  plate.  The  non- 
uniform  inlet  velocity  distribution  was  replaced  by  a  uniform  one. 

The  modified  configuration  representing  a  jet  exhausting  from  a 
plate  at  twice  the  free- stream  velocity  is  shown  in  Fig.  4.  Several 
problems  were  apparent  with  this  initial  configuration.  Addition  of  the 
flat-plate  grid  severely  increased  the  turn-around  time  for  obtaining 
digital  computer  solutions.  The  expanding  nature  of  the  plate  grid  spac¬ 
ing  resulted  in  poor  simulation  of  a  flat  plate,  since  leaks  occurred  in 
some  areas.  With  the  vortex- lattice  program,  pressures  on  a  surface 
can  be  computed  only  at  control  point  locations.  The  flat-plate  grid 
limited  the  pressure  distribution  computations  to  these  points.  One 
possibility  for  improvement  was  to  replace  the  expanding  octagonal 
grid  with  a  much  tighter  rectangular  one  with  uniform  spacing.  How¬ 
ever,  this  solution  would  greatly  increase  the  computer  time  required. 

A  better  solution  involved  using  the  symmetry-plane  capability  of  the 
computer  program.  Since  no  flow  is  allowed  to  pass  through  the  plane 
of  symmetry,  it  provides  an  ideal  infinite  plate  representation. 

A  verification  study  was  conducted  to  prove  the  validity  of  the  sym¬ 
metry  application  to  this  particular  problem.  Figure  5  is  a  schematic 
of  the  three  configurations  which  were  examined.  The  configurations 
were  intended  to  analytically  duplicate  the  experiment  conducted  by 
Mosher  in  the  Georgia  Institute  of  Technology  9-ft-diam  wind  tunnel 
(Ref.  2).  As  Fig.  5  indicates,  the  region  of  interest  included  the  lower 
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surface  of  the  plate  and  the  volume  surrounding  the  jet.  Although  an 
inlet  was  simulated  on  the  side  of  tne  plate  away  from  the  jet.  its  effect 
upon  the  cross-hatched  region  was  deemed  negligible,  assuming  the 
plate  to  be  infinite.  The  model  shown  in  Fig.  5b  was  investigated  to 
define  the  inlet  effect  upon  the  flow  beneath  the  plate  and  to  determine 
the  feasibility  of  shutting  off  the  flow  through  the  vortex  tube.  As  was 
expected,  computations  indicated  that  pressures  on  the  underside  of 
the  plate  were  unchanged  except  near  the  outside  edges.  The  absence 
of  flow  exhausting  from  the  vortex  tube  jet  efflux  simulation  resulted 
in  slight  changes  in  the  flow  field  near  the  tube  exit;  however,  no  var¬ 
iation  was  observed  within  the  region  of  interest.  The  final  configura¬ 
tion,  illustrated  in  Fig.  5c,  produced  results  very  similar  to  those  of 
the  model  shown  in  Figs.  5a  and  b.  Use  of  the  symmetry  plane  to  re¬ 
place  the  vortex  grid  removed  the  plate  edge  effect  entirely.  In  addi¬ 
tion,  removal  of  the  grid  reduced  to  zero  the  disturbances  generated 
by  the  individual  vortex  elements.  Since  no  other  changes  could  be  de¬ 
tected  in  the  region  of  interest,  the  symmetry  plane  model  was  chosen 
as  the  configuration  to  be  used  in  the  subsequent  investigation. 

With  respect  to  the  relative  orientation  of  the  individual  vortex 
panels,  the  model  illustrated  in  Fig.  5c  is  identical  to  the  free  propel¬ 
ler  model  presented  in  Ref.  3.  However,  use  of  the  mirror  image 
plate  representation  requires  omission  of  the  inlet  grid  and  the  asso¬ 
ciated  tube  exit  ring  along  with  a  rotation  of  the  entire  vortex  network. 
This  is  necessary  because  the  computer  program  is  limited  to  provid¬ 
ing  for  symmetry  about  the  Y  =  0  plane.  It  should  be  noted  that  only 
the  input  portion  of  the  vortex  network  is  shown  in  the  figure.  The  re¬ 
flected  portion  of  the  mathematical  representation  is  handled  internally 
by  the  program. 


SECTION  IV 
PARAMETRIC  STUDY 


Upon  completion  of  the  vortex-lattice  model  preparation,  an  at¬ 
tempt  was  made  to  evaluate  the  new  configuration  by  comparing  analyt¬ 
ically  generated  data  with  experimental  results.  In  an  effort  to  more 
closely  simulate  the  experimental  phenomena  and  test  the  sensitivity 
of  the  model  under  various  conditions,  changes  were  made  in  a  number 
of  parameters.  The  parameters  selected  included  width  and  shape  of 
the  trailing  vortex  sheet  and  entrainment  rate  through  the  vortex  tube 
wall.  Information  gained  in  such  an  investigation  would  be  of  value  in 


5 


AEOC-TR-73-57 


the  extension  of  the  jet  representation  to  other  velocity  ratios.  The 
model  used  as  a'basis  for  this  portion  of  the  study  is  shown  in  Fig.  6. 

A  jet- to-free- stream  velocity  ratio  of  eight  was  used  to  define  the  tube 
trajectory.  The  boundary  conditions  controlling  entrainment  through 
the  tube  walls  were  set  at  zero.  The  effective  jet  expansion  rate  and 
width  of  the  trailing  sheet  remained  the  same  as  those  used  originally 
on  the  propeller  configuration.  In  order  to  accomodate  the  change  in 
jet  trajectory,  however,  the  vortex  spacing  along  the  tube  was  modified. 
The  distance  from  the  jet  exit  to  the  first  tube  ring  was  increased  to 
two  and  one-half  times  the  nondimensional  value  used  in  the  propeller 
configuration.  The  parameter  which  controls  the  spacing  increment  of 
the  remaining  tube  ring  stations  was  increased  by  a  factor  of  five. 

These  changes  were  made  to  lengthen  the  lattice  tube  and  remove  the 
tube  exit  from  the  region  of  interest.  In  addition  to  these  and  other 
geometrical  changes  in  the  model,  computations  were  made  for  a  vari¬ 
ety  of  entrainment  rates. 

Five- digit  numbers  were  used  to  identify  the  various  models  in¬ 
vestigated  during  this  study.  The  significance  of  each  digit  is  given  in 
Table  I,  Appendix  II,  along  with  the  implication  of  various  digit  values. 
The  number  "1.  0000"  was  chosen  to  represent  the  model  described  a- 
bove.  Figure  7  illustrates  the  configurations  investigated  along  with 
the  corresponding  model  identification  numbers. 

Three  types  of  analytical  data  were  correlated  with  experimental 
information.  These  were  pressure  distributions  and  streamlines  on 
the  flat  plate  surface,  and  velocity  vectors  at  selected  points  in  the 
flow  field  about  the  jet.  The  results  of  the  investigation  are  described 
in  the  following  paragraphs. 

Analytical  solutions  for  the  initial  configuration,  Model  1.  0000, 
produced  flow  characteristics  considerably  different  from  those  ob¬ 
served  by  Mosher.  Although  a  contrarotating  vortex  pair  was  present 
in  the  analytical  flow  field,  the  sense  was  reversed  from  that  of  the 
vortices  within  a  real  jet.  Addition  of  an  entrainment  velocity  to  the 
tube  surface  boundary  conditions  resulted  in  significant  changes  in  the 
associated  flow  field.  As  is  shown  by  the  vectors  in  Fig.  8,  the  sense 
of  the  vortex  pair  was  corrected.  The  circulation  was  increased  to  the 
extent  that  fluid  was  forced  through  the  tube  lattice.  This  is  illustrated 
by  the  streamlines  on  the  plate  surface  shown  in  Fig.  9  for  three  en¬ 
trainment  values.  Note  that  the  leakage  produces  streamlines  similar 
to  those  generated  by  a  two-dimensional  doublet- sink  comb  inat  ion  jex^ 
cept  in  the  region  downstream  of  the  jet  tube.  The  high  velocities 
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associated  with  the  strong  recirculation  are  reflected  in  the  pressure 
distribution  on  the  plate  surface,  as  can  be  seen  in  Fig.  10. 

On  the  basis  of  this  information,  it  appeared  that  any  improvement 
in  the  analytical  flow  characteristics  would  require  reduction  of  the  re¬ 
circulation  being  generated  by  the  vortex  pair.  Previous  experience 
had  shown  that  this  recirculation  could  be  reduced  by  reducing  the  trail¬ 
ing  vortex  sheet  width.  The  shape  of  the  sheet  attached  to  the  x.  xOxx 
series  models  was  identical  to  that  used  on  the  propeller  jet  configura¬ 
tion  described  in  Ref.  3.  Although  the  propeller  jet  had  been  developed 
for  a  velocity  ratio  of  2.  14,  no  attempt  had  been  made  to  modify  the 
sheet  contour  for  the  plate-jet  combination,  since  no  suitable  expirical 
data  could  be  found  for  a  velocity  ratio  of  8.  0.  In  the  absence  of  such 
information,  an  intuitive  guess  was  made  as  to  the  proper  sheet  con¬ 
tour.  It  is  known  that  rollup  of  the  jet  occurs  farther  along  the  trajec¬ 
tory  for  increasing  velocity  ratios.  To  more  closely  represent  the 
rollup  phenomenon,  it  was  decided  to  elongate  the  vortex  sheet  shape. 

As  shown  in  Fig.  11,  the  high  recirculation  was  reduced.  In  fact, 
the  streamlines  are  seen  to  closely  resemble  those  about  a  two-dimen¬ 
sional  potential  cylinder.  The  plot  of  the  plate  pressure  distribution 
shown  in  Fig.  12  was  also  somewhat  promising.  The  pressures  in  the 
regions  fore  and  aft  of  the  jet  were  much  reduced,  although  the  aft  zero 
pressure  line  was  moved  outward,  further  increasing  the  size  of  the 
positive  pressure  region.  The  extreme  low  pressures  on  either  side 
of  the  jet  were  overrelieved  and  ended  up  somewhat  higher  than  the  ex¬ 
perimental  values. 

Encouraged  by  the  reduction  in  the  recirculation  characteristic  and 
the  favorable  plate  pressure  distribution,  the  authors  decided  to  resume 
the  increase  in  entrainment  rate  through  the  jet-tube  surface.  Plate 
surface  streamlines  generated  for  increased  values  of  entrainment  are 
shown  in  Fig.  13.  The  trend  suggests  that  an  even  higher  value  is  in 
order  if  the  experimental  surface  pattern  is  to  be  matched.  Examin¬ 
ation  of  the  surface  pressure  distribution,  however,  indicates  that  the 
increase  in  entrainment  value  from  Model  1.  2100  and  Model  1.  4100 
was  already  too  large  for  this  particular  configuration.  Figure  14 
shows  that  the  pressures  on  either  side  of  the  tube  are  lower  than  the 
experimental  data.  It  is  interesting  to  note  that  the  increase  in  entrain¬ 
ment  caused  a  universal  decrease  in  pressure  over  the  plate  and  that 
both  positive  pressure  regions  were  drastically  reduced  in  size.  Figure 
15  shows  the  effect  of  further  decreasing  the  entrainment.  It  is  evident 
that  excessive  circulation  is  produced  in  the  jet  cross  section.  This 
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results  in  velocities  considerably  higher  than  those  measured  near  the 
real  jet. 

Model  1.4110  was  generated  in  an  attempt  to  improve  the  plate 
streamlines  without  degrading  either  the  flow  field  below  the  plate  or 
the  plate  pressure  distribution.  As  Figs.  16  and  17  show,  the  attempt 
was  unsuccessful.  The  flow  field  exhibited  the  same  severe  recircula¬ 
tion  and  leakage  through  the  jet-tube  surface  which  had  previously  ap¬ 
peared  in  Model  1.  3000. 

As  a  final  attempt  to  improve  the  analytical  jet  representation  by 
(  reduction  of  the  leakage  through  the  jet-tube  surface,  the  x.  xlOl  series 
of  models  was  generated.  The  new  configuration  was  identical  to  the 
x.  xlOO  model  series  except  for  a  closer  spacing  of  the  vortex  grid  a- 
long  the  jet  tube.  Both  the  initial  space  width  and  the  spacing  incre¬ 
ment  were  reduced  by  50  percent.  Since  no  additional  vortex  filaments 
were  added  to  the  model,  however,  the  modification  resulted  in  a  short¬ 
er  jet  tube. 

The  changes  in  configuration  resulted  in  a  reverse  in  sense  of  the 
circulation  within  the  jet  cross  section.  The  effect  can  be  seen  both 
in  the  plate  streamlines  and  in  the  flow  field  around  the  jet,  as  shown 
in  Figs.  18  and  19,  respectively.  Computations  indicated  that  the  re¬ 
versal  occurred  for  a  wide  range  of  entrainment  rates.  Although  it  has 
not  been  demonstrated,  the  authors  feel  that  the  addition  of  vortices  to 
the  model  was  not  directly  responsible  for  the  problem.  Modification 
for  tighter  vortex  spacing  significantly  shortened  the  tube.  Since  the 
computer  program  automatically  trails  the  vortices  to  infinity  after  the 
last  vortex  ring,  a  severe  change  occurs  in  the  effective  jet  trajectory. 
This  is  shown  schematically  in  Fig.  20.  In  this  configuration  the  tube 
exit  and  the  associated  trailing  vorticity  fall  very  near  the  region  of 
interest.  The  situation  is  further  complicated  by  the  fact  that  the  dis¬ 
continuity  occurs  near  the  point  of  maximum  jet  bending  rate. 


SECTION  V 
CONCLUSIONS 


During  the  course  of  the  investigation  a  number  of  conclusions  were 
drawn  prior  to  the  choice  of  the  basic  vortex  configuration.  Model 
1.  0000.  Many  have  been  discussed  previously  and  are  reiterated  here 
for  the  convenience  of  the  reader.  Others  are  included  for  the  first 
time.  Conclusions  are  as  follows: 
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1.  Preliminary  investigation  indicated  that  the  plane- of- 
symmetry  technique  could  be  successfully  applied  to 
the  jet-plate  combination.  This  was  possible  for  two 
reasons:  first,  the  flow  characteristics  within  the  jet 
tube  were  of  no  interest,  and  second,  an  infinite  flat 
plate  was  being  simulated.  Use  of  this  technique  sim¬ 
plified  the  analytical  modeling,  permitted  a  reduction 
in  computer  time  required,  and  provided  a  means  of 
eliminating  the  complications  involved  in  representing 
the  plate  with  a  vortex  grid.  It  should  be  pointed  out 
that  general  analytical  modeling  situations  involve  in¬ 
lets,  curved  surfaces,  and  multiple  jets.  Consequently, 
the  user  is  deprived  of  the  luxury  of  using  the  symme¬ 
try  plane  to  represent  the  surface  from  which  the  jet 
emerges.  In  many  cases,  however,  symmetry  does 
exist  in  the  flow  field,  and  time  can  be  saved  through 
use  of  the  plane- of- symmetry  technique. 

2.  In  situations  where  inlets  are  present,  caution  should 
be  exercised  in  the  positioning  of  the  vortex  tube  exit. 
Flow  expanding  from  the  tube  exit  can  severely  inter¬ 
fere  with  the  surrounding  flow  field,  especially  when 
high  jet  velocities  are  to  be  simulated.  Therefore, 
the  tube  should  be  extended  sufficiently  far  down¬ 
stream  to  render  the  effects  of  the  expanding  fluid 
negligible. 

3.  Initial  computations  for  Model  1.  0000  revealed  that 
the  trailing  sheet  circulation  was  opposite  in  direction 
from  that  existing  within  a  real  jet.  It  was  shown  that 
the  addition  of  entrainment  through  the  tube  surfaces 
could  correct  the  situation.  The  requirement  for  ad¬ 
ditional  entrainment,  however,  resulted  in  a  severe 
recirculation  and  in  leakage  through  the  tube  walls. 
Similar  conditions  could  be  generated  by  widening  the 
trailing  sheet,  as  shown  by  Model  1.4110. 

4.  With  respect  to  the  flat-plate  pressure  distribution, 

it  was  found  that  high  entrainment  values  could  result 
in  excessively  low  pressures  on  either  side  of  the  jet. 
The  study  indicated  that  a  nearly  identical  pressure 
distribution  was  produced  each  time  the  recirculation 
occurred. 
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5.  The  attempt  to  eliminate  the  recirculation  by  tighten¬ 
ing  the  jet-tube  lattice  spacing  proved  to  be  unsuccess¬ 
ful.  However,  the  reappearance  of  the  reverse  cir¬ 
culation  was  thought  to  be  caused  by  the  shortening  of 
the  tube  rather  than  by  the  closer  spacing  of  the  vor¬ 
tices. 

6.  Of  the  various  configurations  investigated.  Model 

1.  4100  provided  the  best  correlation  with  experimen¬ 
tal  data.  Comparisons  of  the  analytical  and  empiri¬ 
cal  results  are  presented  in  Figs.  21,  22,  and  23. 
Vectors  throughout  the  flow  field  surrounding  the  jet 
efflux  were  in  excellent  agreement.  However,  based 
upon  the  trends  established  by  variation  of  entrain¬ 
ment,  the  optimum  plate  pressure  distribution  would 
occur  at  a  slightly  lower  value  than  that  used  in  Model 
1.  4100.  With  respect  to  streamlines  on  the  plate  sur¬ 
face,  it  appears  that  a  higher  entrainment  value  would 
have  improved  the  data.  This  is  more  obvious  if  the 
streamlines  interpreted  from  the  oil-trace  data  shown 
in  Fig.  5  are  modified  to  remove  the  effect  of  the  pre¬ 
viously  mentioned  contrarotating  eddies  near  the  plate 
surface.  The  modified  streamlines  would  represent 
the  best  possible  results  obtainable  using  an  analytical 
model  which  did  not  include  a  provision  for  simulating 
the  eddies.  If  the  disturbances  did  not,  in  fact,  extend 
far  from  the  plate  surface,  the  streamlines  might  close 
ly  resemble  the  flow  field  immediately  past  the  eddies. 
Figure  24  is  the  authors'  depiction  of  the  proposed 
streamlines. 

7.  The  discrepancies  apparent  in  the  model  1.4100  data 
suggest  that  a  configuration  modification  is  in  order. 

It  is  the  opinion  of  the  authors  that  an  attempt  should 
be  made  to  simulate  both  the  wake  region  aft  of  the 
jet  plume  and  the  pair  of  eddies  near  the  plate.  How¬ 
ever,  no  suitable  technique  is  envisioned  for  a  poten¬ 
tial  flow  model  such  as  the  one  discussed  in  this 
document. 
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SECTION  VI 
RECOMMENDATIONS 


Model  1.  4100  is  capable  of  providing  a  quality  of  data  which  is  far 
superior  to  that  generated  by  most  other  analytical  jet  representations. 
However,  it  is  not  adequate  for  applications  which  require  precise  re¬ 
production  of  the  near  flow  field.  The  following  comments  are  included 
in  order  to  provide  a  starting  point  for  further  development  of  the  plane- 
of- symmetry  technique. 

The  classical  method  of  improving  any  model  generated  using  fi¬ 
nite  element  theory  involves  increasing  the  number  of  elements.  A 
point  of  diminishing  returns  is  reached,  however,  when  the  improve¬ 
ment  realized  becomes  smaller  than  the  errors  associated  with  the 
mismatch  between  the  body  being  represented  and  the  potential  model. 
This  is  not  to  mention  the  drastic  increases  in  computer  time  associ¬ 
ated  with  the  addition  of  vortex  filaments.  It  is  futile  to  exchange  an 
octagonal  vortex  tube  for  one  with  twenty  sides  when  the  latter  is  cap¬ 
able  of  no  more  precise  simulation  of  the  real  flow  characteristics  (i.  e.  , 
viscous  effects,  wake,  etc.  )  than  is  the  former.  The  user  might  find 
it  advantageous,  however,  to  choose  the  octagonal  tube  over  one  of 
square  cross  section  after  examining  the  respective  distances  of  sing¬ 
ularity  error  propagation  into  the  far  flow  field.  It  is  the  opinion  of 
the  authors  that  the  octagonal  tube  used  in  this  investigation  is  adequate 
for  representating  a  real  jet  if  the  region  of  interest  does  not  include 
the  volume  within  one- half  the  diameter  of  the  tube  surface.  Tighter 
spacing  along  the  jet  axis  should  provide  better  simulation  if  additional 
vortices  are  added  to  maintain  the  original  tube  length.  The  authors 
feel  that  the  optimum  model  would  be  a  twelve- sided  tube  with  the  vor¬ 
tices  spaced  evenly  along  the  jet  axis  in  such  a  manner  that  each  indi¬ 
vidual  vortex  panel  is  nearly  square. 

After  the  original  modeling  was  carried  out,  several  improvements 
were  added  to  the  AEDC  vortex-lattice  computer  program.  Since  the 
modified  revision  of  the  program  does  not  limit  "wing  part"  size,  the 
analytical  model  could  be  greatly  simplified  by  substituting  one  pair  of 
"wing  parts"  for  the  two  pairs  currently  used  to  make  up  the  upstream 
half  of  the  jet  tube.  Such  a  modification,  however,  would  require  a 
change  in  the  method  for  distributing  the  "extra  points"  which  define  the 
trailing  vortex  sheet.  A  detailed  description  of  the  method  is  included 
as  Ref.  3. 
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The  trends  established  during  the  course  of  this  investigation  sug¬ 
gest  that  an  even  more  narrow  vortex  sheet  than  that  attached  to  Model 
1.  4100  might  provide  more  realistic  flow-field  data.  Such  a  modifica¬ 
tion  would  probably  necessitate  an  increase  in  entrainment  value,  since 
the  strength  of  the  trailing  sheet  vorticity  would  more  than  likely  be  re¬ 
duced.  It  would  also  be  of  interest  td  define  vortex  sheet  shapes  using 
the  experimental  data  presented  in  both  Refs.  2  and  4  and  to  compare 
them  with  the  contour  associated  with  Model  1.  4100.  In  addition,  it 
might  be  profitable  to  examine  the  flow  about  a  configuration  with  the 
sheet  collapsed  against  the  side  of  the  vortex  tube.  The  authors  have 
investigated  such  a  model  for  the  zero- entrainment  case,  and  it  is  de¬ 
scribed  in  Ref.  3.  The  resulting  flow  field  was  found  to  be  identical  to 
that  produced  by  an  inlet -potential  cylinder  combination,  showing  poor 
correlation  with  experimental  jet  data. 

Since  the  trend  for  a  flat-plate  jet  model  with  entrainment  seems 
to  be  toward  the  more  narrow  trailing  sheet,  i.e.,  closer  to  the  plain 
cylinder  case,  there  is  a  possibility  that  the  real  flow  fields  may  be 
governed  by  different  phenomena.  It  should  be  emphasized  that  in  or¬ 
der  to  maintain  similar  flow  conditions  an  increase  in  entrainment  was 
required  when  the  width  of  the  trailing  sheet  was  reduced.  A  more 
likely  conclusion  is  that  the  two  potential -flow  mechanisms  involved, 
namely  tube  surface  entrainment  and  the  trailing  vortex  pair,  are  re¬ 
sponsible  for  the  flow-field  similarity.  The  addition  of  entrainment  to 
the  propeller  model  described  in  Ref.  3  might  be  a  means  of  providing 
an  insight  into  the  situation. 

Another  avenue  which  has  been  virtually  unexplored  by  the  authors 
involves  the  application  of  an  entrainment  distribution  around  the  jet 
tube,  both  around  the  periphery  and  along  the  axis.  Although  no  suit¬ 
able  empirical  data  seem  to  be  available,  it  may  be  possible  to  deter¬ 
mine  such  distributions  by  computer  experimentation.  Various  distri¬ 
butions  can  be  chosen  and  the  resulting  flow  fields  compared  to  those 
provided  by  experiment. 

Investigation  of  this  approach  should  be  extended  to  include  other 
jet-to-free  stream  velocity  ratios  in  an  attempt  to  define  the  effects  of 
that  variable  upon  the  analytical  configuration.  Such  an  effort,  how¬ 
ever,  would  require  more  empirical  data  than  is  presently  available. 

One  additional  recommendation  is  offered  since  it  is  somewhat  re¬ 
lated  to  the  technique  described  in  this  report.  The  results  previously 
presented  indicate  that  a  potential  model  with  entrainment  can  produce 
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flow  fields  similar  to  those  generated  by  a  configuration  with  trailing 
vorticity.  Jet  cross-sectional  shapes  are  obtainable  from  a  number 
of  existing  experimental  sources  (e.  g. ,  Refs.  2  and  8).  In  addition, 
three-dimensional  analytical  schemes  such  as  the  one  presented  by 
Hackett  and  Miller  in  Ref.  5  are  available.  Based  on  the  present 
study,  a  vortex- lattice  tube  with  its  shape  derived  from  one  of  the 
above-mentioned  sources  and  with  a  suitable  entrainment  distribution 
over  the  surface  could  provide  improved  analytical  simulation  of  a 
real  jet  plume. 
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a.  X/R  =  6.000 

Fig.  1  Jet  in  Crossflow  q-Vector  Projections,  V(/V„  =  8  (Ref.  2) 
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b.  X/R  =  10.000 
Fig.  1  Continued 
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c.  X/R  =  16.000 
Fig.  1  Continued 
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d.  Y/R  =  0 
Fig.  1  Continued 
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e.  Y/R  =  2.000 
Fig.  1  Continued 
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f.  Y/R  =  4.000 
Fig.  1  Continued 
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g.  Y/R  =  6.000 
Fig.  1  Continued 
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i.  Z/R  =  8.000 
Fig.  1  Continued 
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j.  Z/R  =  12.000 
Fig.  1  Continued 
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I.  Z/R  =  20.000 
Fig.  1  Continued 
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b.  Authors'  Interpretation  of  Streamlines 
Fig.  3  Concluded 
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c.  Symmetry  Plane,  No  Flow  through  Tube 
Fig.  5  Jet  Exhausting  from  a  Flat  Plate 
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Fig.  6  Model  l.xOOO 
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Coordinate  system  shoym  In  this  figure  was 
chosen  for  convenience  in  using  computer 
program.  Analytical  data  In  all  other 
figures  Is  presented  in  coordinate  system 
usod  by  Mosher  (Ref.  2). _ 
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a.  X/R  =  6 

Fig.  8  Comparisons  of  Analytical  q-Vector  Projections  from  Models  1.0000  and  1.3000 
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Fig.  9  Comparison  of  Computed  Streamlines  for  Models  1.1 
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Fig.  10  Comparison  of  Computed  Pressure  Coefficient  Contours  for  Models  1.2000  and  1.3000 


AEDC-TR-73-57 


Model 


Fig.  12  Comparison  of  Computed  Pressure  Coefficient  Contours  for  Models  1.2000  and  1.2100 
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Fig.  14  Comparison  of  Computed  Pressure  Coefficient  Contours  for  Models  1.2100  and  1.4100 


AEDC-TR -73-57 


Point 

X/R 

Y/R 

Z/R 

4/ q» 

1.4100 

1.5100 

1 

6 

1 

24 

0.634 

3.389 

2 

6 

6 

24 

1.173 

2,375 

3 

6 

10 

24 

1 .214 

1 .695 

4 

6 

2 

10 

1  .863 

6,136 

5 

6 

10 

12 

1 .740 

2.822 

6 

6 

6 

5 

1.481 

2.250 

7 

6 

10 

5 

1.416 

1.977 

Scale  Model 

q/qM  =  1  - ►  1.4100 

q/qm  -  10 - •>  1.5100 


Fig.  15  Comparison  of  Analytical  q- Vector  Projections  from  Models  1.4100  and  1.5100 
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Fig.  16  Comparison  of  Computed  Streamlines  for  Models  1.4100  and  1,4110 
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a.  X/R  =  6 

Fig.  17  Comparison  of  Analytical  q- Vector  Projections  from  Models  1.4100  and  1.4110 
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Fig.  17  Continued 
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Fig.  19  Continued 


56 


Plane  of  Symmetry 


AEDC-TR -73-57 


Point 

X/R 

Y/R 

Z/R 

q/q«. 

Mosher 

1.4100 

1 

6 

1 

24 

0.95 

0.634 

2 

6 

6 

24 

0.93 

1.173 

3 

6 

10 

24 

1.02 

1.214 

4 

6 

2 

10 

3.51 

1.863 

5 

6 

10 

12 

1.30 

1.740 

6 

6 

6 

5 

1.21 

1.481 

7 

6 

10 

5 

1.20 

1.416 

Scale  Model 


i  »  q/q^,  -  1  - o»  Mosher 


Z 


1.4100 


Intersection  of  Jet 
Boundary  and  X/R  -  6 
Plane  (Ref.  2) 


Vortex  Rotation 


Boundary 


a.  X/R  =  6 

Fig.  21  Comparison  of  Experimental  Data  and  Model  1.4100 
Analytical  q-Vector  Projections 
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Ref.  2  (Mosher) 
Model  1.4100 
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Fig.  24  Proposed  Real  Streamlines  Immediately  Above  the  Surface  Disturbances 
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TABLE  I 

KEY  TO  MODEL  IDENTIFICATION  NUMBERS 


Model 

Number 

First 

Digit, 

Vv" 

Second  Digit, 
Entrainment 
Rate 

Third  Digit, 
Effective  Jet 
Expansion 
Rate 

Fourth  Digit, 
Vortex  Sheet 
Width 

Fifth  Digit, 
Vortex  Ring 
Spacing 

1.  0000 

8.0 

0.  0 

(1) 

(1) 

(2) 

1..1000 
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1. 2000 

1  1 

(1) 

1.  2100 

1.  3000 
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1.  3100 
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1.  3101 
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1.  4100 

0.  64 

(2) 

1. 4110 

0.  64 

(5) 

(2) 

1.5100 

1.  30 

(1) 

(2) 

LEGEND 

(1)  Same  as  contour  used  in  propeller  configuration  (Ref.  3) 

(2)  <  0.  9  jet  diam  from  jet  exit  to  first  ring,  distance  increased  by 

0.  24  jet  diam  for  each  successive  ring 

(3)  Original  contour  elongated  by  2. 1  x  (V-V,,,) 

J 

(4)  0.  45  jet  diam  from  jet  exit  to  first  ring,  distance  increased  by 
0.  12  jet  diam  for  each  successive  ring 

.  (5)  Original  contour  widened  by  200  percent 
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